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SUMMARY 

The optimum sensor location problem, OSLP, may be thought of in terms of the 
set of systems, S, the class of input time functions, I, and the identification 
algorithm (estimator) used, E, Thus, for a given time history of input, the 
technique of determining the OSL requires, in general, the solution of the opti- 
mization and the identification problems simultaneously.* However, this paper 
introduces a technique which uncouples the two problems. This is done by means 
of the concept of an efficient estimator fpr which the covariance of the parameter 
estimates is inversely proportional to the Fisher Information Matrix. 

INTRODUCTION 

The problem of structural identification in structural engineering is one 
which has received considerable attention from several resear ^'ers in the recent 
past (Refs. 1-4). Though various methods have been developed for identifying the 
different parameters that characterize a structure from records obtained in them 
under various loading conditions, few Investigators, if any, have looked at the 
question of where to locate sensors in a structure to lequu. data for "best** 
parametric identification (Ref. 5). The problem of optimally locating sensors in 
a si.ructural system arises from considerations of* (1) minimizing the cost of 
instrumentation; and (2) efficiently detecting structural changes in the system 
with a view to acquiring improved assessment of structural Integrity. 

The problem addressed in this paper can be stated as follows: Given m sen- 

sors, where should they he located in a structure so that records obtained from 
those locations yield the **best’* estimates of the unknown parameters? 

In the past, the optimal sensor location problem (OSLP) was solved by 
positioning the giv<»n number of censors in the system, using the records obtained 
Jit those locations with a specific ectimator, and repeating the procedure for 
different sensor locations. The set of locations which yield the *‘best** parameter 
estimates would then be selected as optimal. The es.:linates obtained, of course, 
would naturally depend upon the type of estimator ujed. Thus the optimal loca- 
tions are estimator dependent, and an exhaustive search needs to be performed for 
each specific estimator. Such a procedure, besides being highly computationa ^ 
intensive, suffers from the major drawback of not yielding any physicrl inc*! 
into why certain locations are preferable to others. 

Recently, work r the solution of the OSIP was done by Shah and Udwadia 
(Ref. 5). In brief, .hey used a linear relationship between small perturbations 
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in a finite dimensional representation of the system parameters and a finite 
sample of observations of the system tLae response. The error in the parameter 
estimates are minimized, yielding the optimal locations. In this paper, we develop 
a more direct approach to the problem which is both con^)ur:ationally superior, and 
throws cinsiderable light on the rationale behind the optimal selection process. 

We uncouple the optimization problem from the identification problem using 
the concept of an efficient estimator (e.g. , the maximum likelihood estimator as 
tine becomes very large). For such an timator the covariance of the parameter 
estima es is a minimum. Using this technique and motivated by heuristic argu- 
ments, a rigorous form lation and solution of the OSLP is presented. The method 
is applied to a building structure modelled as a general linear dynamic system. 
For the N degree of freedom system considc^-ed, the methodology for selecting 
m(m<N) of the nodal displacements for purpostc of measurement is presented. 

Sanq>le calculations are made for a simple building st^cture modelled as a 
two-degree-of-f reedom system subjected to base excitations. The sensor 

location for the identification of: (a) the mass ratio; and (b) the stiffness 

ratio is investigated. 

The results indicate that the OSLP depends on: 

1) the class of systems, S, to which the structure belongs; 

2) the type of excitation; 

3) the actual system parar liters involved; and 

4) the parameters to be identified. 

THEORY 

Consider a systejn modelled by the equation 

MX + CX + KX = F(t) (1) 

where M, C, and K are the (NXN) mass, damping and stiffness matrices, F(t) is 
an (NXl) vector containing inertial forces and extei/^ally applied loads and X 
is the N-vector of nodal displacements. Let 0j^, 6(3 and be vectors containing 
^he various parameters related to the mass matrix, the, danq)ing matrix and the 
stiffness matrix, respectively, which need to be Identified. For convenience, 
we collect these quantities in the parameter vector, 0, defined as 



where the superscript T indicate matrix transpose. If the M, C and K are 
symmetric each of the three subvectors has a maximum dimension of N(N+l)/7. 

Given m sensors (m < N) , we jfhen need to find where to locate them so that 
the covariance of the estimate, o, is a minimum. Assulu^ further that the 
measurement vector Z(t) can be expressed as 
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i = 1,2, ...,N 


( 2 ) 


= g, [X(e,t)] -h N^(t) 

where Zj^ is the ith component of Z(t), and the functionals gj^ represent the 
"measurement process*’- The dependence of the response X on the parameter 
vector 0 is explicitly noted The measurement noise N^(t) is taken as non- 
stationary Gaussian White noise with a variance of i|/2(t). Therefore, 

E[N\(t^)N^(t2)]= . (3) 

where 6j^ and 6j) stand for the kroneker and the dirac-delta functions, 
respectively. A total of m out of N responses need to be selected so that 
they contain the most information about the system parameters and are maxi- 
mally sensitive to any changes in the parameter values. This "selection" 
process can be represented by an m-dimensional vector Y such that 

Y(t) = 3Z(t) (4) 

where S is the (m ^ N) upper triangular selection matrix with each row con- 
taining null elements except for one which is unity. The m different compo- 
nents of Z selected to be measured are so ordered in vector Y, that if the 
element in the i-th row and k-th column of S is unity, the (1+1) -1th row has 
unity in its Jt-th column with t > k. The matrix S has the property that 
P = s^s in an (NXN) diagonal matrix with unity in its i-th row if, and only 
if, Zj^ is selected to be measured* The elements of P are otherwise zero. 

Hence, one can write 

Y(t) = Sg[X(9,t)] + SN(t) (5A) 

4 HfX(0,t)] + V(t) (5B) 


If is linearly related to che response X^, in general, then 


H(X(6,t)] = SRX 


( 6 ) 


where R(t) can be thought of as a dynamic gain matrix* In the case that 
is relat^’ to the resporse X^ only, then matrix R will reduce to a diagonal 
matrix, o . J. 

The problem of locating sensors in an optimal manner then reduces to 
determining the selection matrix S, or alternatively, finding the m locations 
in P chat 'should be unity. These locations must be so chosen as to obtain the 
"best" parameter estimates. 

SOME MOTIVATING THOUGHTS AND THE FISHER INFORMATION MATRIX 


Consider a »^ase in which one tries to estimate only one parameter, 0x (io 
be identified) involved in a dynamic system model with only one sensor provided. 
Therefore, one wants to ideally choose a location i (out of N possible such 
locations) such that the measurement y±(t)^ ie[l,N], te(0,T) at location i yields 
the best estimate of the parameter 0i. Heuristically, one should place the 
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sensor at such a location that the time history of measurements obtained at that 
location is most sensitive to any changes in the parameter 9]^. He^^ce, in equa- 
tion (5B) it is really the slope of H(X(0i,t)J with respect to 8i that needs to 
be maximized. However, since only the absolute magnitude of this slope is of 
interest, it is logical to want t find i (or equivalently determine the selec- 
tion matrix S described pteviously) such as to maximize (3H/36x)^ over the inter- 
val (0,T) during which the response is to be measured. This leads to maximizing 
the following integral: 



When there is more than one parameter to be estimated, and the number of 
sensors is greater than unity, this intuitive approach needs to be extended in a 
more rigorous iianner. In such cases recourse to mathematical treatment is 
necessary, and we shall see that such treatment will be in agreement with our 
heuristic solution outlined above. 

To further understand the problem, let us loo'c at it from another angle, 
namely, the concept of an efficient unbiased estimator. For such as estimator 
che covariance of the estimates is a min imum . Furthermore, it can be shown that 
for any unbiased estimator of 9, 


(9-0)(9-0)'^ > 

/(Hf 

(II) 

1 J 

Lo 



( 8 ) 


where § is the estimate of 0 and the matrix [9H/30]x^ 4 3Hx/30j. If the esti- 
mator is *’ef f icient", the above inequality becomes an equality. This means that 
the left-hand side of inequality (8) takes its lowest value (minimum covariance) . 
Hence, 


E (0-0) (0-0)" 



(9) 


The term inside the bracket on the right-hand side of the equation (9) is known 
as the Fisher Information Matrix, Q(T). Thus, maximi’^ing Q(T) would indeed lead 
to a minimization of the covariance of the estimate, 0. 

We note then that the m sensor locations need to be so chosen that a suitable 
nortr l the matrix Q(T) given by 


Q(T) 



0 


( 10 ) 
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is maximized. This constitutes an extension of equation (7), which we heuris- 
tically derived earlier for the scalar case, to the vector situation. Introducing 
equation (6) in equation (10) one may write 


Q(T) 


/ 


T T 


dt. 


( 11 ) 


where the ij element of can be written as: 

L J j 

where X = and 0 = note that the Fisher Matrix is symmetric and 

is dependent on che length of the record available, as well as the locations of 
the sensors as determined by the matrix P. 

If the m locations where the sensors are to be placed are denoted by 
Sj^, k = l,2,,.-,ni, then 


P = 


m 


E 

k=l 



( 12 ) 


where the (N x N) diagonal matrix* Is^ has all its elements equal to zero except 
the element of the row, which is unity. Noting that P is a diagonal matrix, 
equation (11) can be simplified to yield 


Q[T;3^,s^,...,s^;S,0;I] = 


m ^ 

E / 

k-l . 


T T 

T X„dt 




(13) 


where is the row of the matrix R. Also in eq. (13) explicit mention Is 
made of the dependence of the Fisher !iatrlx on the time length T of the available 
data, the syscem S, the parameter vector 6, and the time-variant input I. If the 
matrix R is diagonal, with diagonal elements then the ij element of 

the matrix Q, after some manipulation, reduces to 


^ -f 


9x 3 
s. s, 
k k 

30 ^ 36 ^ 



(14) 


Each element of represents the cross-sensitivity of measurement with respect 

to the response x of node s, . 

s t ^ 

k 
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The cptimal sensor locations are then obtained by picking m locations sj^, 
k = l,2,,..,ni, out of a possible N, so that a suitable norm of the matrix Q is 
maximized (e.g., the trace norm, etc...)- This may be specified by the 
condition 

max I I Q |T 5 Sn , » • • • 1 s 5 S , 0 j [ { • (15) 

e(l,N) L ® J 

Although there are several matrix norms which could be used, perhaps the 

most useful and physically meaningful in this context is the trace norm. In 
order not to detract the reader from the basic methodology we defer an exhaustive 
treatment of suitable matrix norms to a later communication. 

The methodology presented up to this point is valid for both linear and non- 
linear systems since the criterion developed in equation (13) was derived using 

only equations (5) and (9). We will now indicate its application to linear 
multi-degree-of-f reedom systems. 

APPLICATION TO LINEAR DYNAMIC SYSTEMS 

Consider the N-degree-of -freedom dynamic system whose governing differential 
equation of motion is given by eq. (1), together with X(tQ) = Xq> X(tQ) = 
where Xq and Xq are the given initial conditions for the system. Assume the 
system to be classically damped. Introducing 


X(t) = <^n(t) 

where ♦ is the (N x N) weighted modal matrix and ri(t) is the N-vector of 
generalized coordinates we get 

o 

"n + = dp(t), n(tQ> » n (t^) « MX^, 

where the (N x N) diagonal matrix A is given by 


M • ■ [''xl- 

The solutic.i of equation (17) is given as 

t 

hi(t) = Hq u^(t-t )+fiQ v^(t-tQ)+ f h^(t-T)p^(T)dx 

i i J 

0 

where and are initial conditions and 


(13) 


(17) 


(18 
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Coscjj t + Sinoj, t 

'•i '*1 


Uj<t) . EXp(-{.Uj^t) 

”l‘"> ■ TiT Eip(-?iV>)=‘"“d.'' 




hi(t) = v^(t). 


“d. " “i\A - . and 


p.(c) = i = 1,2,. ...N. 

Also, differentiating equation (1) with respect to 0, yields 

/ /\ \ a 

MXg+CXg+KXg = Fg(t) - (MQX+CgX+KQXj; Xq(0) = 0, Xq(0) = 0 


where 


K]« - 5 ^ 


with 


A. • 

M X = [M X : M X : M X ; M X] 

e • °2 j L 


X drid j * 

Introducing 


Xe = 


yields 

z+25jjOij^z+Az = G(t) 

where 


G(t) - «’^[f,-(m,X+C^X«,x)J 


Equation (21) can further be simplified to give 

■/\ A A 


G(t) - <t> |^FQ-(Mg$n+CQ$rr+Kg$njJ 


(19) 


(20) 


( 21 ) 


( 22 ) 
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where f| and n can be obtained by differentiation of eq- (18). This may be shown 
as follows 


n^(c) = y* '5^(t-T)pj(T)dT 


where 


W^(t; = -EXPj 


Y^(t) = EXp(-5^a.^tj 



(A) . + 

1 1 i / 

■*1 “dj 


Sinti) L , 


Cosu)j t 


\\) 


H^(t) = and 

P^(t) = $“F(t). i = 1,2,. 


,,N. 


AliO 


t 

n^(t) = / 


hj(t-T)pj^(x)dT 


where 



(23) 


(24) 
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Y^(t) 







Sino)j t 


- 25 ^(d.Cosu)j t 
i 1 

h,(t) = Y,(t), 

X i 

P^(t) = $'^F(t), i=l,2,...,N. 



Therefore, substituting equations (23) and (24) into equation (22) gives G(t). 
Consequently the solution of equation (20) can be written as: 


z..(t) = f h (t-t)G. . (i)dT 
ij J i ij 

0 


(25) 


where hi(t) is the same as that of eq. (18). Notice that the initial conditions 
in eq. (20) are zero. This is due to the fact that the initial conditions of 
(18) are known constants. 

If we assume that [C] is expressed as a linear combination of [K] and [M] , 
then eq. (22) can further be simplified. Namely, 

C = 2aK+2BM, (26) 


where a and 3 are known constants. Hence in equation (17), the percentage of 
damping, can be expressed as: 


^1 = 


oicj. + 

1 0)1 


' y i 1 , 2 ^ , N 


(27) 


To further simplify equation (22) under this assumption, let us consider the 
following three cases: 

1) The vector 0 contains only i.e., only estimation of mass param- 
eters is undertaken. Then 

T • 

G(t) ^ ^ [F^ - M^iJ)(n+26Ti)] . (28A) 

2) The vector 0 contains only the subvector 0jr. Then 
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1 . 

G(t) = 4 . [Fg - Kg(Kn+2an)] 


(28B) 


3) Finally if the vector 0 = [a g] , 


X • T • X 

G(t) “ F - 2 An, (p F - 2In^ 

QL B 


(29) 


Lf the input F(t) is not a function of 6 , than p 0 would be omitted all through 
this discussion. Once the solution of equation (25) is obtained, the Fisher 
iiatrices may be obtained as in equation (13). Hence 


T T T 

T z r r 4>z 
m /• S- s, 

Q - dt 


E/ 

k=l •' 

^ 0 




(30) 


Je note that the summation form of relation (30) is particularly amenable to the 
maximization of the trace norm of Q. 
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EXAMPLE 


To illustrate some of the ideas of the previous section, consider the problem 
of finding the optimal sensor location (OSL) in a structural system modelled by the 
two-degree-of -freedom system (shown in Figure 1) which is subjected to the base 
excitation of f(t). 

The governing differential equation of motion can oe expressed as 


M + CX + KX = -W f(t) 


(31) 


T T 

where X = <x^ x^> , C = aK, W = <Am and the matrices M and K are 



- 

- 


- 



A 

0 


B+1 

-1 

M = 

0 

1 

m , and K = 

-1 

1 


A case study for locating sensors to best identify (1) the mass ratio, A, of the 
first to the second floor and (2) the stiffness ratio, B, of the first to the 
second floor, will be presented. 

Let s][ denote the lower mass location and S2 the upper mass location. The 
selection between the locations can be equated to determining the one non-zero 
element of the [1x2] selection matrix, S, with the measurement H(t) defined by 

H(t) = SX + V(t) , 


where, V^t) is Stationary Gaussian White Noise (S G W N) with i(/(t)=i/^^. 


If S = [1 0] the lower mass is selected for measurement; if S « [C 1] the 
upper mass is selected. The location si would then be preferred over the loca- 
tion S2 for identifying the parameter A, if Q[T,S}^] > Q[T,S2], where T is the 
time that the measurement is taken. 



(32 A) 
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Since only one parameter is being estimated the Fisher matrices reduce to 
scalars . 

'fhe dependence of the OSL on various types of the base excitations can be 
studied now. Let us for this presentation consider ground acceleration in the 
form of a delta function, i.e., f(t) = <5(t) . 

In this case, closed form solutions foi and can be obtained. 

For the OSL problem for the **best" (minimum covariance) identification of 
the parameter A (given the parameters B and a) using an Impulsive base excitation. 
Figure 2-A shows the plots of tha ratio of the information matrices Ql(T)/Q 2 (T), 
for T = 50 secs, for various values of the parameters A (which is to be identi- 
fied) and a* 4 o^^o> where u)o 4 Points on the graph with ordinates greater 

than unity indicate the optimal location to be the lower mass level and vice 
versa. The graphs indicate that the optimal location In most cases, for the 
range of A considered, is the upper mass level. However, we observe that for 
some small values of A and a* the OSL is the lower level. We note, interestingly 
enough, that the optimal sensor location for identification of A actually depends 
not only on the actual values of B and a which are presumably known, but also on 
the value of the parameter A itself which is to be identified! Thus to be able 
to ascertain the optimal sensor location some a priori assessment of A is 
necessary. 

Figure 2~B shows that the optimal location for identification of the param- 
eter B (given A and a), using an Impulsive base input, is again the upper mass 
level for the range of B values considered. For larger B values, however, and 
a*>0,05, the trend appears to be more and more in favor of the upper mass. This 
seems intuitively correct, for as B becomes larger, the lower part of the system 
becomes stiffer and the OSL would be the upper mass level. 

Figure 2-C is associated with the OSLP for estlrating the parameter B usir : 
a sinuscii^al base excitation, f(t) =* a sin cjt. The figure shows that as the 
normalized driving frequency y = w/wo varies, the OSL changes. For this example 
the Fisher Matrices can be computed in closed form. For the e stimat ion of B, 
(g''ven A and a* =* 0) the dimens ii.nl.ess driving frequency y = /1+1/A yields no 
information on B from records at either of the two mass levels. The responses at 
the two mass levels yield identical amounts of information on B at y = 0 and 
y = /l for k ^ 1, as indicated by the values of Q1/Q2 - 1 at these frequencies. 
The value of Qi/Q 2 = 0 at y = 1 is indicative of the fact that the upper mass 
level is . far better location for a senior when estimating B with a* = 0. 

Figure 2-D shows the mean value of the r«.tio Q1/Q2 a random Gaussian white 
noise base excitation together with the 1-a band. The OSL appears to be at the 
upper mass level for identification of A. 
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CONCLUSIONS AND DISCUSSION 


This paper presents a general raethodolog> for determining the optimal sensor 
locations in dynamic systems for obtaining records which would enable the **best** 
(minimum covariance) identification of a given set of unknown parameters in the 
system. The technique utilizes the concept of an efficient estimator to uncouple 
the identification from the optimization problem. In order to present the basic 
idea in as clear a fashion as possible, we have restricted the discussion In this 
sequel tc linear systems. 

The method has been illustrated by application to a two degree of freedom 
system. Though the results presented here for the simple system chosen form only 
a first step towards acquiring a detailed understanding of the OSL problem, the 
following conclusions appear to be relevant at this time: 

(1) The OSL for a given system heavily depends on the class of forcing 
functions used for obtaining response data. In this study, an im- 
pulsive base motion Is considered. 

(2) The OSL for linear dynamic systems is Independent of the amplitude of 
the forcing function. 

(3) The OSL depends in general on all the values of system parameters. For 
instance, the OSL for estimating A with minimum covariance depends not 
only on the actual parameter values B a but on the value of A it- 
self for the system! This implies that the OSL problem associated with 
identifying a given parameter (or a set of parameters) in a dynamic 
system necessitates the knowledge of some a priori estimates of the 
unknown parameter (s) . 

(4) ihe results of our simple example show that che OSL problem may yield 
solutions which may be diificult to predict on pureJy heuristic 
grounds. The OSL appears to depend, even fur this relatively simple 
problem, in a rather complex manner on the actual parameter values of 
the system and the nature of the base excitation. 
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Figure 2-A* Variation of Q 1 /Q 2 va.^ous 
values of the parameter A. 
Q 1 /Q 2 greater than unity 
Indicates that optimal is 
at lr*w*ir mass. 


Figure 1. 2^degree--of-f reedom 
gener'c structural 

syste.-n. 



Figure 2-B. Variation of Q 1 /Q 2 for various 
values of B. 
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